Upper limit on the photon fraction in highest-energy cosmic rays from AGASA data 
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A new method to derive an upper limit on photon primaries from small data sets of air showers 
is developed which accounts for shower properties varying with the primary energy and arrival 
direction. Applying this method to the highest-energy showers recorded by the AGASA experiment, 
an upper limit on the photon fraction of 51% (67%) at a confidence level of 90% (95%) for primary 
energies above 1.25 ■ 10 20 eV is set. This new limit on the photon fraction above the GZK cutoff 
energy constrains the Z-burst model of the origin of highest-energy cosmic rays. 
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PACS numbers: 96.40.Pq,96.40.-z,13.85.-t,13.85.Tp 

Since their first discovery about 40 years ago Q, the 
existence of particles with energies around and above 
100 EeV = 10 20 eV was confirmed by several air shower 
experiments using different measurement techniques 2, 
03- tJ S 0- Q| ■ The quest for the nature and origin of 
these extremely high-energy (EHE) cosmic rays contin- 
ues to drive considerable experimental and theoretical 
efforts B. As first pointed out by Greisen, Zatsepin and 
Kuzmin |lj], the travel distance of EHE particles is lim- 
ited due to energy losses on background radiation fields. 
For instance, the energy loss length of 200 EeV protons is 
~30 Mpc (e.g. ISII3)- A cosmological origin of the ob- 
served EHE particles thus seems disfavoured. However, 
no astrophysical object in our cosmological vicinity could 
be identified yet as source of these events. Moreover, ex- 
plaining efficient particle acceleration to such enormous 
energies poses a theoretical challenge. The acceleration 
problem is circumvented if EHE particles are generated 
in decays or annihilation of topological defects (TD) or 
super-heavy dark matter (SHDM) [lj [II El [II . These 
objects are expected in certain inflation scenarios and 
have also been proposed as dark matter candidates. 

A common feature of such non-acceleration models is 
the large fraction of EHE photons in the injected parti- 
cle spectrum. Due to interactions of these photons with 
background fields, the diffuse photon flux at GeV en- 
ergies allows one to derive an upper limit on the elec- 
tromagnetic energy injected as EHE particles at dis- 
tances beyond a few Mpc [TH fl6| . This constrains non- 
acceleration models which predict particle injection at 
large distances |l2l Il7j . In turn, models with injection 
sites closer to the observer imply a significant fraction of 
primary photons in the observed EHE events. As an ex- 
ample, in the SHDM model metastable particles of mass 
M x ~ 10 14 GeV are clumped in the galactic halo and pro- 
duce EHE photons, nucleons and neutrinos by decay |12fl . 
Thus, stringent limits on the EHE photon fraction pro- 
vide constraints on non-acceleration models complemen- 
tary to those from the GeV photon background. 

Based on an analysis of muons in air showers observed 
by the Akeno Giant Air Shower Array (AGASA), up- 
per limits on the photon fraction were estimated to be 
28% above 10 EeV and 67% above 32 EeV (95% CL) [|. 



Comparing rates of near- vertical showers to inclined ones, 
upper limits of 48% above 10 EeV and 50% above 40 EeV 
(95% CL) were deduced from Haverah Park data Q. 
Non-acceleration models of cosmic-ray origin are not 
severely constrained by these bounds 18], however. Pho- 
tons are predicted to reach a considerable fraction only 
at highest energies, while with decreasing energies below 
100 EeV the "conventional" hadronic cosmic-ray compo- 
nent soon outnumbers photon primaries due to the steep 
flux spectrum. For instance, based on the SHDM model 
the photon fraction above 40 EeV is estimated as ~25% 
only, increasing to ~50% at 70 EeV 12]. 

In this work, we focus on events above 100 EeV. These 
particles are most directly linked to the production sce- 
nario in non-acceleration models, and the predicted pho- 
ton dominance can be checked with the data. The largest 
data set on EHE events available to date was obtained 
by the AGASA experiment. From eleven AGASA show- 
ers reconstructed with energies above 100 EeV, the muon 
content in the shower is measured in six 0, ITsj . For each 
event, adopting its reconstructed primary parameters, we 
compare the observed muon signal to results from shower 
simulations for primary photons. In contrast to the anal- 
ysis method in 5], where the data distribution above 
energies of 10 EeV and 32 EeV was compared to an over- 
all simulated distribution, we thus use the information 
about the individual event characteristics. We develop 
a new statistical method that allows us to combine the 
information from all events and to set a limit on the pri- 
mary photon contribution. 

AGASA 0,113] consisted of 111 array detectors spread 
over ~100 km 2 area and 27 muon detectors with an en- 
ergy threshold of 0.5 GeV for vertically incident muons. 
The primary energy was determined with a statistical 
accuracy of ~25% for hadron primaries |l9j. Assum- 
ing photon primaries, the energies reconstructed this 
way were found to be underestimated by ~20% for the 
most-energetic events Six events were reconstructed 
with >100 EeV which had more than one muon detector 
within 800-1600 m distance from the shower core @. The 
muon density pj at 1000 m core distance was obtained 
for each event j=1...6 by fitting an empirical lateral 
distribution function plj to the data. The uncertainty 
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TABLE I: Reconstructed shower parameters of the AGASA 
events [f| (upper part of the Table) and simulation results 
(lower part). The energies are increased by 20% to account 
for the case of photon primaries The azimuth angle is 
given clockwise from north for the incoming direction. 



primary energy [EeV] 


295 


240 


173 


161 


126 


125 


zenith angle [°] 


37 


23 


14 


35 


33 


37 


azimuth angle [°] 


260 


236 


211 


55 


108 


279 


Pj [m" 2 ] 


8.9 


10.7 


8.7 


5.9 


12.6 


9.3 


preshower occurrence [%] 


100 


100 


96 


100 


93 


100 


<P S j> [m -2 ] 


4.2 


3.1 


2.1 


2.3 


1.7 


1.8 


1.1 


1.0 


0.9 


0.6 


0.5 


0.5 


x 2 


1.6 


3.0 


3.4 


2.2 


4.6 


4.0 


20.8 


8.3 


6.4 


13.9 


3.1 


4.6 



estimated for the resulting pj is 40% The recon- 
structed shower parameters of the highest-energy events 
with muon data are given in Tab. Q] 

It is well known that photon-initiated showers gen- 
erally contain significantly fewer muons than hadron- 
induced events. For each AGASA event, 100 primary 
photon showers were generated. The reconstructed pri- 
mary direction |19| was chosen as simulation input and 
the primary energy varied from shower to shower accord- 
ing to the reconstruction uncertainty. The energy was 
also globally increased by 20% to account for the energy 
underestimation in case of photons. Electromagnetic cas- 
cading of photons in the geomagnetic field was simu- 
lated for the AGASA site with the new PRESHOWER 
code 22] . For the AGASA events, in most cases 
preshower formation occurred (see Tab. [IJ. The re- 
sulting at mosp heric shower was simulated with COR- 
SIKA 6.18 [23] as a superposition of subshowers initiated 
by the preshower particles or, if no preshower occurred, 
with the original primary photon. Electromagnetic inter- 
actions were treated by the EGS4 code |24|, which was 
upgraded [23] to take photonuclear reactions as well as 
the Landau-Pomeranchuk-Migdal (LPM) effect |2f| into 
account. For the photonuclear cross-section, we chose the 
extrapolation recommended by the Particle Data Group 
(cr PDG ) [HI23 shown in Fig.[I] The influence of assum- 
ing different extrapolations is discussed below. Hadron 
interactions were simulated with QGSJET 01 [U, and 
for energies < 80 GeV with GHEISH A [13. 

The distribution p* of simulated muon densities ob- 
tained from CORSIKA for each AGASA event is given 
in Fig. |2 together with the data. The average values 
<pj> and standard deviations Ap S j are listed in Tab. [Q 
The average muon densities for primary photons are a 
factor 2-7 below the data. To quantify the level of agree- 
ment between data and primary photon expectation, a 
X? value is calculated for each event j as 
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with the measurement uncertainty Apj = 0.4 ■ pj [5j. To 
account for possible deviations of the simulated muon 
densities from a Gaussian distribution, the probability 
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FIG. 1: Data |2(! and extrapolations of the photonuclear 
cross-section a Jp . The PDG extrapolation (er PDG ) |2r3 . l27j| 
is chosen for this analysis. Also shown are two parametriza- 
tions with larger cross-sections at highest energies, denoted 

^mod ^extr gj| ^ rpj^ cross _ section on air ig 
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FIG. 2: Observed muon densities (points with error bars) 
compared to the muon densities expected for primary photons 
(histograms) for the six events. Assigned to each event the 
primary energy (see Tab. HJ. The measured muon densities 
are larger than predicted for primary photons. 



Pj(.X 2 ^Xj) °f a photon-initiated shower to yield a value 
X 2 > X 2 is determined by a Monte Carlo technique: A 
simulated muon density is taken at random from the dis- 
tribution p s j , a random shift is performed according to the 
experimental resolution Apj , and a % 2 value is calculated 
with Eq. (Q, replacing pj by the artificial muon density 
value. Repeating this many times gives ft (x 2 >x 2 ). The 
values xj and Pj are listed for the six events in Tab. [I] 
The probabilities pj range from 3% to 21%. 

The combined probability £>(x 2 >X^=i xf) °f s ^ x 
photon-initiated events to yield a x 2 value larger or equal 
to the measured one is p — 0.5%. Thus, it is unlikely that 
all cosmic rays at these energies are photons (rejection 
with 99.5% confidence), and it is possible to derive an 
upper limit on the primary photon fraction F^. 

It should be noted that, due to the small event statis- 
tics, the upper limit cannot be smaller than a certain 
value. Assuming a primary photon fraction F y , a set of 
n m primaries picked at random ab initio does not contain 
any photon with probability (1 — F^) nm . For n m =6, this 
probability is ~5% for F 7 =40%. Thus, in the present 
case only hypothetical photon fractions F 7 >40% could 
in principle be tested at a confidence level of 95%. 

For deriving an upper limit F" 1 <100%, scenarios have 
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FIG. 3: Upper limits (95% CL) on cosmic-ray photon frac- 
tion derived in the present analysis (P) and previously from 
AGASA (A) ,5] and Haverah Park (H) 3 compared to some 
predictions from super-heavy dark matter (SHDM) [T^ . Z- 
burst (ZB) and topological defect (TD) |l3fl models. 

events might be initiated by photons. For a hypothet- 
ical photon fraction F 7 , the probability q that a set 
of n m events contains n 7 photons is q{F 1 ,n^ ll n ra ) — 



F" 7 (l - F 7 ) n = 



,(n B 



This probability is multiplied 



by the probabilities p 7 (n 7 ) • pj(n m — n 7 ), with p 7 (n 7 ) 
being the probability that the n 7 most photon-like look- 
ing events are generated by photons, and Py(n m — n 7 ) 
being the probability that the remaining n m — n 7 events 
are due to non-photon primaries. p 7 (rt 7 ) is determined 
by the MC technique as the probability to obtain val- 
ues x 2 > xl t i with p 7 (0)=l and with x\=X 2 from 
Tab.[I] where the index k\ refers to the event with small- 
est value Xj i an d xli — X% 1 ■ To derive an upper limit on 
photons, the probabilities py(n m — n 7 ) are set to unity. 
Summing over all possibilities n 7 =0 . . . n m then gives the 
probability -P(-F 7 ) to obtain x 2 values at least as large as 
found in the data set, 



P(F-y) = z] g(F 7 ,n 7 ,n m ) ■ p 1 {n 1 ) ■ p^ (n B 



(2) 



This probability depends on the assumed photon frac- 
tion F 1 . For the considered data set one obtains 
P(F 7 =51%) = 10% and P(F 7 =67%) = 5%. There- 
fore, the upper limit on the primary photon fraction is 
F 7 ' = 51% (67%) at 90% (95%) confidence level. 

The derived bound is the first experimental limit on the 
photon contribution above the GZK cutoff energy. The 
limit refers to the photon fraction integrated above the 
primary photon energy that corresponds to the lowest- 
energy event in the data sample, which in the present 
analysis is 125 EeV. In Fig. [21 upper limits derived 
previously at lower energy and the current bound are 
compared to some predictions based on non-acceleration 
models. Models predicting photon dominance at highest 
energies are disfavoured by the new upper limit. 

The statistical stability of the upper limit can be tested 
by, e.g., omitting one event and calculating an upper limit 
with the remaining five. Iterating through all six possi- 
bilities of rejecting an event, the upper limits are between 
61-78% (95% CL). Alternatively, the 320 EeV Fly's Eye 
event can be added to the event list with the photon 



probability of 13% 33]. The upper limit then is 66% 
(95% CL). The result is quite stable, as the individual 
photon probabilities do not differ much from each other. 

The upper limit derived in the present analysis is con- 
servative with respect to different sources of systematic 
uncertainties, since pj might be overestimated. As an 
example, in the 295 EeV event data, muon detectors sat- 
urated and the obtained pj might rather be regarded 
as a lower limit pj, 0|. Concerning the simulations, 
p S j is robust when changing the low-energy hadronic in- 
teraction model |35| . The applied high-ener gy model 
QGSJET 01 produces -20-30% more muons ; 36j com- 
pared to SIBYLL 2.1 and also compared to a prelim- 
inary version of QGSJET II j^. Smaller values of p s j or, 
in case of the 295 EeV event, a larger value of pj, would 
decrease pj and reduce the photon upper limit. 

The derived upper limit is robust against reasonable 
variations of the primary photon energy adopted. In gen- 
eral, a larger primary energy would result in larger values 
of pj and pj. We already accounted for a possible 20% 
underestimation in case of primary photons. It seems un- 
likely that an additional, systematic underestimation of 
reconstructed primary photon energies of more than 20- 
30% exists, also because of the stronger preshower effect 
at increased energy that makes the profiles of primary 
photon showers more similar to hadron-initiated events. 
In turn, a rescaling of AGASA energies to smaller values 
would make the muon densities predicted for photons 
even more discrepant to the data. 

A considerable uncertainty exists in extrapolating the 
photonuclear cross-section. A stronger (weaker) increase 
of the cross-section with energy than adopted in this 
work leads to larger (smaller) values of p s j. We repeated 
the calculations with different cross-section assumptions. 
The upper limit of 67% (95% CL) changes little for mod- 
est variations of the extrapolation. Adopting, for in- 
stance, the parametrization denoted cr mod in Fig. ^ the 
upper limit becomes 75% (95% CL). However, as an illus- 
tration, assuming the extreme extrapolation labeled <7 extr 
(Fig. , the simulated p s j are increased on average by 70- 
80% with respect to the calculation using er PDG . In such 
a scenario, one would obtain P(cr cxtr , F 7 =100%)~ 15%, 
and no upper limit could be set with high confidence. 
Also theprevious limits from Haverah Park and AGASA 
data 0,O| would increase when assuming cr oxtr . 

In summary, we introduced a new method for deriving 
an upper limit on the cosmic-ray photon fraction from 
air shower observations. Applied to the highest-energy 
AGASA events, an upper limit of 67% (95% CL) is ob- 
tained for cosmic rays >125 EeV. This photon bound im- 
poses constraints on non-acceleration models of cosmic- 
ray origin, with possible implications also on the descrip- 
tion of the dark matter or inflation scenarios in these 
models. Within the next few years, a considerable in- 
crease of EHE event statistics is expected from the HiRcs 
detector [|| and the Pierre Auger Observatory (3i|. Thus, 
even more stringent conclusions on EHE photons are pos- 
sible by performing an analysis as presented here. It will 
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be studied elsewhere, to what extend a scenario of dom- 
inant EHE photons together with a large photonuclear 
cross-section can be tested with shower data. 
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